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ABSTRACT 

In a previous letter, we computed the decay constants of heavy pseudoscalar mesons in the 
framework of relativistic (instantaneous) Bethe-Salpeter method (full 0“ Salpeter equation), 
in this letter, we solve the full 1“ Salpeter equation and compute the leptonic decay constants 
of heavy-heavy and heavy-light vector mesons. The theoretical estimate of mass spectra of 
these heavy-heavy and heavy-light vector mesons are also presented. Our results for the 
decay constants and mass spectra include the complete relativistic contributions. We find 
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1 Introduction 


The decay constants of mesons are very important quantities [1, 2, 3, 4] and the study of them has become 
an interesting topic in recent years [5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18]. These constants play 
important roles in many aspects, such as in the determination of Cabibbo-Kobayashi-Maskawa matrix 
elements, in the leptonic or nonleptonic weak decays of mesons, in the neutral D — D or B — B mixing 
process, etc. 

In a previous letter [19], the decay constants of heavy-heavy and heavy-light pseudoscalar mesons are 
calculated in the framework of relativistic instantaneous Bethe-Salpeter method [20] (also called Salpeter 
method [21]), good agreement of our predictions with recent lattice, QCD sum rule, other relativistic 
model calculations as well as available experimental data is found. 

In this letter, we extend our previous analysis to include vector mesons, present the relativistic 
calculation of heavy-heavy and heavy-light vector decay constants in the framework of full Salpeter 
equation. The instantaneous Bethe-Salpeter equation, which is also called full Salpeter equation, is a 
relativistic equation describing a bound state. Since this method has a very solid basis in quantum field 
theory, it is very good in describing a bound state which is a relativistic system. In a earlier paper [22], 
we solved the full 0“ Salpeter equation of pseudoscalar mesons; in another paper [23], we solved the full 
1 Salpeter equation of equal-mass vector mesons. The predictions of these relativistic methods agree 
very well with other theoretical calculation and recent experimental data. In this letter, we extend the 
analysis of equal-mass vector system to non-equal mass system by solving the full 1“ Salpeter equation 
of vector mesons, and use this method to predict the values of decay constants for heavy-heavy and 
heavy-light vector mesons. There are some input parameters in our method, we need to hx them when 
we solve the full 1“ Salpeter equation. In our calculation, we fix the input parameters by fitting the 
experimental data of mass spectra, so we also present the mass spectra for heavy vector mesons. 

This letter is organized as following, in section 2, we introduce the relativistic Bethe-Salpeter equation 
and Salpeter equation. In section 3, we give the formula of relativistic wave function and decay constants 
of vector meson. We solve the full Salpeter equation, obtain the mass spectra and wave function of vector 
mesons. Finally, we use these relativistic wave function to calculate the decay constants of heavy vector 
mesons and show the numerical results and conclusion in section 4. 


1 



2 Instantaneous Bethe-Salpeter Equation 


In this section, we briefly review the instantaneous Bethe-Salpeter equation and introduce our notations, 
interested reader can find the details in Ref. [22]. 

The Bethe-Salpeter (BS) equation is read as [20]; 

f d'^k 

- mi)xiq){l^2 + m 2 ) = i J ■j^^V{P,k,q)x{k) , (1) 

where xiq) is the BS wave function with the total momentum P and relative momentum q of the bound 
state, and V{P, k, q) is the kernel between the quarks in the bound state. pi,P 2 are the momenta of the 
constituent quark 1 (heavy or light) and anti-quark 2 (heavy or light), respectively. The total momentum 
P and the relative momentum q are defined as; 


Pi = aiP + q, ai = 
P 2 = a 2 P -q, 02 = 


mi 


mi -|- m2 

m 2 


mi -|- m 2 

The BS wave function xiq) satisfy the following normalization condition; 


d^kd'^q 


Tr 


^(Pi)'S '2 ^iP 2 )S'^ik - q) + ViP,k,q)] xiq) 


= 2iPo 


( 2 ) 


where Siipi) and 52 (^ 2 ) are the propagators of the two constituents. In many applications, the kernel 
of BS equation is “instantaneous”, i.e., in the center mass frame of the concerned bound state (P= 0), 
the kernel R(R, k, q) takes the simple form; 


ViP,k,q)^Vik,q) = Vi\ k\,\ 9 \,cos9) , 


where 0 is the angle between the vectors k and q. Then the BS equation reduced to the Salpeter 
equation. 

It is convenient to divide the relative momentum q into two parts, gy and g_L, a parallel part and an 
orthogonal one to the total momentum of the bound state, respectively. 


q^ = q'{\+q^^, 

gf = (P•g/M2)P^ = 

Correspondingly, we have two Lorentz invariant variables; 
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= T > 9t = • 

In the center of mass frame P= 0, they turn out to the usual component go and | 9 |, respectively. Now 
the volume element of the relative momentum k can be written in an invariant form: 


d^k = dkpk^dk^dsdcj) , 


where cf) is the azimuthal angle, s = {kpqp — k ■ q)/{k^qj,). The instantaneous interaction kernel can be 


rewritten as; 


V{\k- a \) = V{k^,qus) . 


Let us introduce the notations (pp{q^ and for three dimensional wave function as follows: 


* J ^xiqpq±), 


d{q±) = J ^^^^^V{k^,q±,s)ipp{k'^) . 


Then the BS equation can be rewritten as: 


x{q\\,q±) = Si{pi)ri{q^)S2{p2) ■ 


The propagators of the two constituents can be decomposed as: 


g.(p.) =_ ^ipiq±) _^_ \p{q±) _ 

J{i)qp + aiM - uii + ie J{i)qp + aiM + tOi - ie ' 


, - IIP 

^i = \Jm‘f + ql, Aj(g^) = ^ ± ^ (8) 

where i = 1,2 for quark and anti-quark, respectively, and J(i) = (—Here A^(g_L) satisfy the 
relations: 


A+(g±) + Kpiq±) = ^ > ^tpiq±)^^tpiq±) = Aj(g±), ^tpiq^)^^Uq±) = o ■ (9) 

Due to these equations, A^ may be considered as P—projection operators, since in the rest frame ^ = 0 
they turn to the energy projection operator. 

Introducing the notations p^^{q±) as: 


v^p^iq±) = ^fp{q±)^Mq^)^^fpiq±) 
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and taking into account with = 1, we have 


( 9 ±) + + v>p {q±) 


With contour integration over Qp on both sides of Eq.(6), we obtain; 


v’p{q±) 


Ajp(g±)??p(g±)A^p(g±) _ ^ipiq±)Vp{Q±)^2pi<l±) 

{M-(jJi-U}2) {M + U)i+0J2) 


and we may decompose it further into four equations as follows: 


(M - wi - + (g±) = A^p(<?_L)^/p(g±)A^p(g_L) , 

(M + wi + u;2)</?p"(g±) = -Kp{qi-)vp{qi-)^2p{qi-) > 

</^p“(g±) = = 0 . (11) 

The complete normalization condition (keep all of the four components appearing in Eq.(ll)) for BS 
equation turns out to be: 


qldqT 

2tx‘^ 


Tr 


^++^(^++— - Tp — ^(p — — 

^ M ^ M 


= 2Pn 


( 12 ) 


3 Relativistic Wave Function and Decay Constant of Vector Meson 


The general form for the relativistic wave function of vector state = 1“ can be written as 16 terms 
constructed by P, q, e and gamma matrix. Because of the approximation of instantaneous, the 8 terms 
with P ■ q± become zero, so the general form for the relativistic Salpeter wave function of vector state 
= 1“ can be written as [23]; 


q^i-{q±) = • ei 


fi{q±) + -^f 2 {q±) + ^f3{q±) + ^^fi{q±) 


+ Af^x/5(<?±) 


+f'±FMq±) + {4 lA - qi. ■ A)Mq±) + ^{FA4 i. - Fqi. ■ A)h{q±)^ 


(13) 


where the is the polarization vector of the vector meson. The equations 


(/9+_ (g±) = ifAiqA = 0 


(14) 


give the constraints on the components of the wave function, so we have the relations 


/i(9±) 


[q±f 3 {q±) + Ar^/ 5 (g±)] (mim2 - U>1CU2 + q'j) , . 

M(mi + m2)qj_ ’ 


f5(q±)M(-mim2 + u;iu ;2 + qj) 
(mi - m2)qj^ 
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f / N [-Qlh{Q±) + M‘^h{q±)] {miuj2 - m 2 U)i) f6{q±)M{miuj2 - ^2^1) 

j 2 {q±) — T\jr,. I ,. ^^2 ’ js{q±) 


M{uJi + uJ2)q2_ (ivi - ui2)qj_ 

Then there are only four independent wave functions f 3 {q±), f 4 :{q±), f 5 {q±) and f6{q±) been left in the 
Eq.(13). Following the Ref.[22], put Eq.(13) into Eq.(ll) and take trace, we obtain four coupled integral 
equations, by solving them we obtain the numerical results of mass spectra and wave functions for the 
corresponding bound states. 

In our calculation, we choose the center-of-mass system of the heavy meson, so gy and g_L turn out 
to be the usual components {qq, 0) and (0, uji = (mf + and 102 = {ni^ + ■ Wave functions 

/ 3 (C)) / 4 (C)) /sC^) and f%{q) will fulfill the normalization condition: 


dq lQuJiU )2 inf f 


+ 


ujii 02 - mim 2 + q 




74/5 - /s + -^6 


= 2M. (15) 


J (27r)3 3 ' (mi + m2)(wi + ^ 2 ) 

In our model, Cornell potential, a linear scalar interaction plus a vector interaction is chosen as the 
instantaneous interaction kernel V [22]: 


A 


Vsiq) = -i- + Vo) 5 Yq) + 


Viq) = Vs{q) + %^^^VYq) , 

A 1 


a 


r2 ,^2 


^ (9 + 

where the coupling constant a* (9) is running: 


, V,{q) = - 


2 as{q) 

371^ 21 

[q +a^) 


as{q) = 


127r 


1 


27 


log(a + 


QCD 


and the constants A, a, a, Vq and Aqcd are the parameters that characterize the potential. 
The decay constant Fy of vector meson is defined as 


(16) 


(0|gi7/i92|E,e) = FvMe^l, 


which can be written in the language of the Salpeter wave functions as: 


(0|gi7M^2|E,e) = VNc I Tr 

Therefore, we have 

Fv = WN, 




(27r) 




dq 

(^ 




(17) 


(18) 


(19) 
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Table 1: Mass spectra in unit of MeV for heavy vector meson. ‘Ex’ means the data from experiments 
[24], ‘Th’ means the predictions from our theoretical estimate. 



B* 

Bt 


B* 

D*s 

n* 

■^U 

Ex(lS) 


5416.6 

5325.0 

5325.0 

2112.1 

2010.0 

2006.7 

Th(lS) 

6336.9 

5416.6 

5326.2 

5322.9 

2112.0 

2010.2 

2006.5 

Th(2S) 

6918.5 

5957.6 

5842.3 

5837.7 

2673.0 

2545.9 

2540.8 


4 Numerical Results and Conclusion 

In our method, there are some parameters that have to be fixed when performing the calculations. In 
Ref. [22], we fixed the values of the input parameters by fitting the mass spectra for heavy pseudoscalar 
mesons of 0“ states, and we hope to choose the same parameters for pseudoscalar and vector mesons. But 
we hnd, when we solve the full Salpeter equation Eq.(ll) of heavy vector mesons with same parameter 
set as in Ref. [22] for pseudoscalar mesons, our predictions of mass spectra of vector mesons do not agree 
well with experimental data, and can not explain the vector-pseudoscalar mass splitting. We argue that 
we choose a very simple interaction kernel Eq.(16), while the forms of pseudoscalar and vector wave 
functions are very different (see Eq.(13) in this letter and Eq.(20) in Ref. [22]). The latter decrease the 
connection between the pseudoscalar and vector states, so we like to choose different parameters to ht 
experimental data and to present the mass splitting between the vector and pseudoscalar mesons. We 
find the following parameters can fit data very well, 

a = e = 2.7183 , a = 0.06 GeV, Eq = -0.49 GeV, A = 0.21 GeV^ Aqcd = 0.27 GeV and 

mb = 5.158 GeV, me = 1.7551 GeV, m* = 0.535 GeV,md = 0.377 GeV, mu = 0.371 GeV, (20) 

and we show our predictions of mass spectra for heavy vector mesons as well as the experimental data 
in Table 1 and Table 2. 

In table 1, we show the results for ground state (IS) and first radial excitation state (2S), one can see 
that our mass predictions for ground states of heavy-light mesons can fit the experimental data [24] very 
well. In table 2, we show the mass spectra of the hrst eight states for vector cc system. As can be seen 
from the table, our mass results below 2D state agree well with experimental data, while the masses for 
2D and 45 states are about 50 ~ 60 MeV lower than the experimental data. 
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Table 2: Mass spectra in unit of MeV for cc vector system. 



IS 

2S 

ID 

3S 

2D 

4S 

3D 

5S 

Ex(cc) 

3096.916 

3686.093 

3770.0 

4040 

4159 

4415 



Th(cc) 

3096.8 

3690.9 

3759.8 

4065.2 

4108.2 

4344.2 

4371.6 

4567.2 


Table 3: Mass spectra in unit of MeV for bb vector system. 



IS 

2S 

ID 

3S 

2D 

4S 

3D 

5S 

Ex(66) 

9460.30 

10023.26 


10355.2 


10580.0 


10865 

Th{bb) 

9460.3 

10029 

10130 

10379 

10438 

10648 

10690 

10868 


We also calculate the mass spectra for vector bb system, we find our prediction with upper parameter 
set Eq.(20) can not fit experimental data. The reason is due to that in bb system, there are double heavy 
b quarks, and the flavor iVj = 4, so we have to choose a new set of parameters as well as smaller value 
of coupling constant [25]. We change the previous scale parameters to A.qcd = 0.21 GeV, ruh = 5.1242 
GeV, and other parameters are not changed. With this set of parameters, the coupling constant at the 
scale of bottom quark mass is as{mb) = 0.23, and obtained the mass spectra of bb systems. We show 
the numerical results and experimental data in Table 3. One can see that our predictions can fit the 
experimental data very well, even for higher states. 

By fitting the mass spectra of heavy mesons, we fixed the parameters and obtained the relativistic 
wave functions for heavy mesons. Put the obtained wave functions into Eq.(19), we calculated the decay 
constants for heavy-heavy and heavy-light vector mesons. In Table 4, we show our estimates of decay 
constants for heavy-light ground state (IS) and first radial excitation state (2S) as well as the B* vector 
mesons. In Table 5, we show our estimates of decay constants for cc and bb systems. We also show the 
theoretical uncertainties of our results for decay constants in Table 4 and Table 5. These uncertainties 
are obtained by varying all the input parameters simultaneously within ±10% of the central values, and 
taking the largest variation of the decay constant. 

Eor comparison, in Table 6, we show our predictions for decay constants and recent theoretical 
predictions as obtained by other methods. Eor example, we show the results from Ref. [30], which is 
also in BS method, but they used different interaction kernel, different form of wave function, especially. 
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Table 4; Decay constants of heavy vector meson in unit of MeV. 



Fb* 

FBt 

Fb* 

Fb* 

Fp. 

Fp* 

d 


IS 

418±24 

272±20 

239±18 

238±18 

375±24 

341±23 

339±22 

2S 

331±21 

246±13 

222±15 

221±14 

312±17 

290±16 

289±16 


Table 5: Decay constants in unit of MeV for cc and bb vector systems. 



IS 

2S 

ID 

3S 

2D 

4S 

3D 

5S 

Fy{cc) 

459±28 

364±24 

243±17 

319±22 

157 ±11 

288±18 

174±12 

265±16 

Fvibb) 

498±20 

366±27 

261±21 

304±27 

155±11 

259±22 

176±10 

228±16 


different reduction method from ours, they chose Thompson equation to reduce the full BS equation, 
while we choose instantaneous approach to reduce the full BS equation. We also show the ratios of decay 
constant Fb*/Fb* and Fbi*/Fb* in Table 6. Not like the pseudoscalar case, where we find good agreement 
for pseudoscalar decay constants between different models, from Table 6, we find rough agreement between 


the values of vector decay constants estimated by different methods, this 

means we need 

more effort for 

the knowledge of vector decay constants. 





Table 6; Decay constants and ratios of decay constants estimated by different methods. (NRQM: Nonrelativistic 
Constituent Quark Model, RQM: Relativistic Quark Model, QL: Quenched Lattice QCD, BS: Bethe-Salpeter 
method, RM: Relativistic Mock Meson Model) 

Ref. Fb* 

-^.s 

Fb* 

Fb* / Fb* 

Fp, 

Fp* 

Fp% / Fp* 

ours 272±20 

238±18 

1.14T0.08 

375±24 

340±22 

1.10T0.06 

NRQM [26] 236tn 

i5i;ji 

1 c:c+0.07 
-‘-•^^-0.06 

326li7 

223li^ 

1 41+0-06 
-‘■•^-‘--0.05 

RQM[27] 214 

195 

1.1 

335 

315 

1.06 

QL[28] 217 

190 

1 .10(2)t2 

254 

234 

1.04(1)12 

QL[29] 229±20lf^ 

196±24;^^ 

1.17(4)1° 

272±16l^o 

245±20l^ 

1.11(3) 

BS[30] 

164 


242 

237 

1.02 

RM[31] 225±9 

194±8 

1.16T0.09 

298±11 

262±10 

1.14T0.09 


There are other interesting quantities, such as the ratios of vector to pseudoscalar decay constant 
Fy/Fp, which are sensitive to the difference between the vector and pseudoscalar wave functions. In 
Table 7 we show our estimates of these ratios. 

In conclusion, we calculated the decay constants of heavy vector mesons in the framework of the rela- 


















































Table 7; Ratios of vector to pseudoscalar decay constant. 


Ft 

Fb* 

Fb* 

Fj/-qi 

^n* 

Fo* 

^Vb 

Fbs 

Fb 


FDs 

Fd 

1.37±0.18 1.30±0.24 

1.26±0.28 

1.21±0.27 

1.57±0.23 

1.51±0.26 

1.48±0.26 


tivistic Salpeter method. Our relativistic estimate results are Fd* k. 375 ± 24 , Fq* ^ 340 ± 23 {D*^, D*^), 
Fsi ~ 272 ± 20 , Fb* ~ 238 ± 18 (5*°, B*^), Fb* ~ 418 ± 24 , Fj/^ ^ 459 ± 28 , F^(2S) ^ 364 ± 24 , 
Fy ~ 498 ± 20 and Fx( 25 ) ~ 366 ± 27 MeV. 

This work was supported by the National Natural Science Foundation of China (NSFC). 
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